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Abstract. We address the problem of uniquely reconstructing F-type icosahedral qua- 
sicrystals from few images produced by quantitative high resolution transmission electron 
microscopy and explain recent results in the discrete tomography of these sets. 



1. Introduction 

Motivated by the request of materials science for the unique reconstruction of quasicrys- 
talline structures from a small number of images produced by quantitative high resolution 
transmission electron microscopy (HRTEM), the discrete tomography (DT) of F-type icosahe- 
dral model sets is concerned with the inverse problem of uniquely reconstructing finite patches 
of these aperiodic point sets from their (discrete parallel) X-rays in several directions. As 
for the discrete analogue of parallel X-rays in computerized tomography (CT), the X-ray 
of a finite subset of Euclidean 3-space in a certain direction assigns to each line parallel to 
this direction the number of points of the set on this line. In fact, for some crystals the 
technique QUANTITEM (Quantitative ANalysis of The Information from Transmission 
Electron Microscopy) can approximately measure the number of atoms lying on lines paral- 
lel to directions that guarantee HRTEM images of high resolution, i.e., yield dense lines in 
the crystal; cf. |X5|, [T7] , Therefore, the classical setting of DT deals with the corresponding 
inverse problem for (periodic) lattices; cf. [12], [9l [10], [TT] . Since F-type icosahedral model sets 
are commonly regarded as good mathematical models for many icosahedral quasicrystals [8] 
and since it is reasonable to expect that future developments in technology will extend the 
technique QUANTITEM to classes of quasicrystals, we feel that it is about time to explain 
recent results in the DT of these sets in a manner that is easily accessible for practitioners. 
In particular, special emphasis is put on illustrating the problems that arise in the DT of 
aperiodic model sets. For detailed proofs and related results, we refer the interested reader 

to guana mi. 

2. F-TYPE ICOSAHEDRAL MODEL SETS 

We denote the golden ratio by r, i.e., r = (1 + v5)/2. Note that r is a root of X 2 — X — l£ 
Z[X], whence it is an algebraic integer of degree 2 over (Q. The unique non-trivial Galois 
automorphism of the real quadratic number field <Q(r) = Q(\/5) = (Q © Qr, determined 
by \/5 i— ► — \/5, will be denoted by .', hence r' = 1 — r. Throughout this text, elements 
of Euclidean 3-space will be written as row vectors. The standard face-centred icosahedral 
module of quasicrystallography is given by 

M F := Z[r](2,0,0)ffiZ[r](r + l,T,l)©Z[r](0,0,2), 
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Figure 1. A few slices orthogonal to (r, 0, 1) of a patch of the icosahedral 
model set A from Example 12.21 (left) and their .*-images inside the icosahedral 
window s + W (right), both seen from the positive x-axis. 

where Z[t] = 7L © Zr is the ring of integers in (Q(r); cf. [2] [7]. Clearly, A*f is a free Z[r]- 
module of rank 3 that spans all of 1R 3 . In particular, it is a free Z-module of rank 6. Moreover, 
A^f has full icosahedral symmetry, i.e., it is invariant under the action of the rotation group 
Y of the regular icosahedron centred at the origin G R 3 with orientation such that each 
coordinate axis passes through the mid-point of an edge, thus coinciding with 2-fold axes of 
the icosahedron. Due to the connection with the icosian ring, we prefer to use the scaled 
version L := ^Aif C Q(t) 3 instead of Aip itself; cf. [7J and references therein. 

Definition 2.1. Let the map .* : L — > R 3 be defined by applying the Galois conjugation .' 
coordinatewise. Given a subset W C R 3 with ^ W° C W C W° and W° compact, and 
any t G R 3 , we obtain an F-type icosahedral model set A(t, W) by setting 

A(t,W) :=t + {ae L\a* eW} 

The map .* : L — > R 3 is the so-called star map of A(t, W) and W is referred to as the window 
of A(t, W). The model set A{t, W) is called generic if it satisfies 914^ n L* = 0. Moreover, it 
is called regular if the boundary dW has Lebesgue measure in R 3 . For a window W C R 3 , 
we denote by Xj(W) the set of generic F-type icosahedral model sets with a window of the 
form s + W, where s G R 3 . 

We refer the reader to [H] for the corresponding cut and project scheme and to [161 E] 
for related general settings and background. F-type icosahedral model sets A are aperiodic 
Delone subsets of 3-space. Moreover, if A is regular, then A is pure point diffractive. If A is 
both generic and regular, and, if a suitable translate of its window has the full icosahedral 
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Figure 2. The central slice S of the patch of A from Figure [U (left) and its 
/-image S* inside the (marked) decagon (s + W)n(r', 0, l) 1 - (right), both seen 
from perpendicular viewpoints. 

symmetry of L*, then A has full icosahedral symmetry in the sense of symmetries of Ll-classes, 
meaning that a discrete structure has a certain symmetry if the original and the transformed 
structure are locally indistinguishable; cf. [31 [T3] and references therein for details. 

Example 2.2. For a generic regular F-type icosahedral model set with full icosahedral sym- 
metry, consider A := A(0, s + W), where s := 10 _3 (1, 1, 1) and W is the regular icosahedron 
centred at the origin G H 3 with orientation such that (r',0, 1) and (—r',0, 1) belong to its 
vertices; see Figured] for an illustration. 

3. Complexity 

The fundamental notion in DT is the following. 

Definition 3.1. Let F be a finite subset of R 3 , let u 6 § 2 be a direction, and let C u be 
the set of lines in direction u in 1R 3 . The {discrete parallel) X-ray of F in direction u is the 
function X U F :£„->N :=NU {0}, defined by 

X U F{£) :=card(FnO = J]l F (x). 

Moreover, the support (X U F)~ 1 (N) of X U F, i.e., the (finite) set of lines in C u which pass 
through at least one point of F, is denoted by swpp(X u F). Further, we denote by £^ the 
subset of C u consisting of all lines in C u which pass through at least one point of L. 

Let U C S 2 be a finite set of pairwise non-parallel directions and let F be a, finite subset of 
R 3 . Clearly, F is contained in its grid with respect to the X-rays in the directions ofU given 

G£:=n( U 

u&U £Gsupp(X u F) 
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Figure 3. The grid G that is generated by the slice S from Figure with 
respect to two L^^'^-directions with slopes and tan(27r/5), respectively. 

In fact, if there are more directions in U than elements in F, one even has F = G^; cf. (141 
Proposition 5.3]. Since after about 3 to 5 images taken by HRTEM, typical (quasi)crystalline 
probes may be damaged or even destroyed by the radiation energy, this result is meaningless in 
practice. Consequently, one has to deal with grids that are much bigger than their generating 
finite sets; see Figure for an illustration in the plane. In fact, for a finite subset F of a 
fixed F-type icosahedral model set, its grid can contain points of a different translate of L 
than F itself; see [H Figure 5] for illustrations of this phenomenon in the case of planar sets. 
As mentioned in the introduction, only directions that yield dense lines in a fixed F-type 
icosahedral model set A are reasonable in view of applications in practice. Therefore, we may 
restrict ourselves to A- directions, i.e., directions parallel to non-zero interpoint vectors of A. 
It turns out that the set of Tl-directions is exactly the set of L- directions (defined analogously) 
for any F-type icosahedral model set A; cf. [14(, Proposition 3.20]. 

Definition 3.2 (Reconstruction Problem). Let W C R 3 be a window and let u±, . . . , u m £ § 2 
be m > 2 pairwise non-parallel L-directions. The corresponding reconstruction problem is 
defined as follows. 

Reconstruction. 

Given functions p Uj : C Uj — > No, j G {1, . . . , m}, whose supports are finite and satisfy 
supp(p Uj ) C , decide whether there exists a finite subset F of an element of lJ(W) that 
satisfies X Uj F = p Uj , j G {1, . . . , m}, and, if so, construct one such F. 

In the case of planar lattices and two lattice directions, the above reconstruction problem 
clearly is intimately related with the problem of reconstructing 0-1-matrices from their row 
and column sums. Though, for aperiodic model sets the problem is much more involved due 
to the complication with the window. Below, an L-direction will be called an L^ 0,1 ) -direction 
if it lies in the hyperplane orthogonal to (r, 0, 1), the latter representing a 5-fold axis of the 
icosahedral symmetry of L. Observing first that generic F-type icosahedral model sets can 
be sliced orthogonal to (r, 0, 1) into (planar) cyclotomic model sets that are based on the 
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Figure 4. The image G* of the grid G from Figure [3] under the star map 
together with S* and the corresponding decagonal window already shown in 
Figure [2 Note that the star map naturally extends to Q(t) 3 which contains 
the grid G. Although the star map may change the directions and the relative 
position of parallel lines, the marginal sums are being preserved. In fact, G* 
is the grid generated by 5* with respect to the directions that have slopes 
and tan(27r/10), respectively. The reconstruction task is to find a subset of 
G* that lies in a translate of L*, fits the marginal sums and lies in a suitable 
translate of the interior of the window. 

ring of integers in the 5th cyclotomic field ([14| Proposition 3.16]), the following tractability 
result can be easily proved by using similar arguments and the same methods from algebra, 
graph theory and computational geometry as in [3], where algorithmic complexities in DT of 
cyclotomic model sets were studied; see Figures [I]H] for illustrations. 

Theorem 3.3. [144 Theorem 4.3] When restricted to two L^' ' 1 ^ -directions and polyhedral 
windows, the problem RECONSTRUCTION can be solved in polynomial time in the real RAM- 
model of computation. 

Remark 3.4. For a detailed analysis of the complexities of the reconstruction problem in 
the case of B-type icosahedral model sets, we refer the reader to [131 Chapter 3]. Note that 
even in the case of planar lattices and the Turing machine as the model of computation the 
corresponding reconstruction problem is NF-hard for three or more lattice directions (defined 
analogous to Tl-directions); cf. [10] • Therefore, it seems to be rather obvious that one cannot 
expect a generalization of Theorem 13.31 to the case of three or more L-directions. 



4. Uniqueness 

In general, the above reconstruction problem can possess rather different solutions. There- 
fore, one is led to the investigation of the uniqueness problem of finding a small number of 
suitably prescribed L-directions that eliminate these non-uniqueness phenomena. 

Definition 4.1. Let £ be a collection of finite subsets of 1R 3 and let U C S 2 be a finite set 
of directions. We say that £ is determined by the X-rays in the directions of U if different 
elements of £ cannot have the same X-rays in the directions of U. 

Below, a finite subset C of a Delone set A C R 3 is called a convex subset of A if it satisfies 
the equation C = conv(C) D A; cf. [9] for the lattice case. The set of all convex subsets of 
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A will be denoted by C{A). Using the slicing of F-type icosahedral model sets into certain 
cyclotomic model sets again, one obtains the following fundamental result. 

Theorem 4.2. [H Theorem 5.12] 

(a) There is a set U C S 2 of four pairwise non-parallel L^ 1- ' ' 1 ) -directions such that, for 
all generic F-type icosahedral model sets A, the set C(A) is determined by the X-rays 
in the directions of U. 

(b) For all generic F-type icosahedral model sets A and all sets U C § 2 of three or less 
pairwise non-parallel lA' ' 1 ) -directions, the setC(A) is not determined by the X-rays 
in the directions of U . 

In fact, for a generic F-type icosahedral model set A, the set C(A) is determined by the 
X-rays in the directions of any set U of four pairwise non-parallel L( r ' 0,1 )-directions with the 
property that there is no U -polygon in A; cf. [13]. Here, for a set U of L( r,0,1 )-directions, a 
[/-polygon in A is a planar non-degenerate convex polygon P with all its vertices in A such 
that any line in ft 3 that is parallel to a direction of U and passes through a vertex of P 
also meets another vertex of P. It turns out that one can choose four pairwise non-parallel 
L( r,0,1 )-directions which provide uniqueness and yield dense lines in F-type icosahedral model 
sets; cf. |14[ Example 5.13 and Remark 5.14] for examples and details. Moreover, there 
are even four pairwise non-parallel L^' ' 1 ) -directions which provide uniqueness, yield dense 
lines in F-type icosahedral model and have the property that, in an approximative sense, 
for any fixed window W C IR 3 whose boundary dW has Lebesgue measure in R 3 , the set 
{J AeI F( W -)C(A) is determined by the corresponding X-rays; cf. [H] for details. This setting, 
which we already used in the definition of the reconstruction problem, is highly relevant 
in practice of quantitative HRTEM because it reflects the fact that, due to the icosahedral 
symmetry of genuine F-type icosahedral quasicrystals, the determination of the rotational 
orientation of a quasicrystalline probe in an electron microscope can rather easily be achieved 
in the diffraction mode. Though, in general, the X-ray images taken in the high-resolution 
mode do not allow us to locate the examined sets. Therefore, as already explained in [41 114j. 
in order to prove practically relevant and rigorous results, one has to deal with the whole 
local indistinguishability class of a regular, generic F-type icosahedral model set A, rather 
than dealing with a single fixed one. 

5. Outlook 

Although the presented results answer some of the basic problems of DT of F-type icosa- 
hedral model sets, there is still a lot to do to create a tool that is as satisfactory for the 
application in materials science as is CT in its medical or other applications. Foremost, this 
is due to the fact that there is always some noise involved when physical measurements are 
taken, whereas the results in this text can only be applied when exact data is given. Therefore, 
it is necessary to study stability and instability results in DT of F-type icosahedral model 
sets in the future pQ. 
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